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TECHNICAL MEMORANDUM 


ELECTRICAL TORQUES ON THE ELECTROSTATIC GYRO 
IN THE GYRO RELATIVITY EXPERIMENT 

I. INTRODUCTION 


The ultimate accuracy of an electrostatic gyro is determined by the 
level to which the Newtonian electrical torques can be reduced. An 
understanding of these torques is necessary for extrapolation of per- 
formance in a one-g environment to that expected in a nearly zero-g 
environment. The success of the Stanford Gyro Experiment depends on 
the reduction of all Newtonian gyro drifts to levels well below the 
relativity drifts, which are expected to be 7 arc sec /yr and 0.05 arc sec/yr 
for the geodetic and motional effects, respectively. It is important, first, 
to determine whether this is feasible and, second, to predict the accu- 
racy to which the relativity effects can be measured, since this is 
important in assessing the relative desirability of performing the gyro 
experiment vis a vis other experimental tests of relativity theory. In 
addition, it is important to know the level of gyro performance on the 
ground at which a successful space experiment can be undertaken. For 
all these reasons, we give a general derivation of expressions for the 
torque on an electrostatic gyro from first principles. We also give a 
complete discussion of the subject , including numerical estimates of the 
torques for current state-of-the-art rotors. The various averaging 
effects are also discussed. 

This report relies heavily on work done by Honeywell [1]. In 
Section II , we sketch the derivation of the basic torque equation derived 
in Reference 1 . In Section III , wc describe the sample rotors used for 
the purpose of our calculation. Section IV gives expressions for the 
torque in terms of the lower harmonics. These expressions are similar 
to those of Reference 1 but are those appropriate for circular electrodes. 
Section V gives the basic method of exact numerical integration we have 
developed and represents the main extension of the method of Reference 
1. Our method includes the effect of all the higher harmonics up to the 
20th. These harmonics are shown to be important and should be included 
in the calculation because some of the symmetry properties of the lower 
harmonics may not persist in the higher ones. Also, the Fourier series 
used in Reference 1 does not converge very rapidly. Numerical results 
for the three sample rotors are given in Section VI. Section VII gives 
results for the special case of a gim baled gyro. This is not the actual 
configuration planned for the gyro experiments, but it illustrates the 
large reduction in torques for the gimbaled case. Section VIII describes 
how the method of Section V can be extended to secondary torques, i.e., 
those due to nonuniformity of the rotor electrode gap. The secondary 
torques are not calculated numerically in this report, but their magnitude 


is estimated. The remaining sections describe the three averaging 
effects: the orbital averaging (Section IX), the averaging due to space- 
craft roll (Section X), and the averaging due to rotor spin (Section XI). 
The orbital averaging combined with spacecraft roll is shown to be 
effective for two of the throe components of gravity gradient forej . and 
the third component cun be minimized by proper design of the spacecraft. 
The averaging due to spacecraft roll is shown to be effective for the 
primary torques and is crucial for the achievement of the goal of 1 
milUarcsec/yr drift for the gyro. The averaging due to spin is described 
but not calculated, since a complete map of the rotor is not yet available. 
Section XII states the conclusions of tho report. 
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II. GENERAL DISCUSSION OF ELECTROSTATIC TORQUES 

For a conductor in an external electric field £ , the force per unit 
area f is given by [1] 



- 12 

where c q = 8. 85 x 10 farads/m. The electric field is perpendicular 

to the surface of the conductor and zero inside the conductor. If we 
consider an irregularly shaped conductor, then the total force F and 
torque T are given by 


*- c iG (i) 

s 


f = Y |2 1 2 (r x n) ds , ( 2 ) 

* s 


where the integral is taken pver the entire surface S, n is a unit vector 
normal to the surface, and r is a vector frqjn some origin to the surface 
element dS. It is important to realize that T is dependent on our choice 
of origin. The evaluation of electric torques on a gyro is accomplished 
by choosing a suitable model for the actual shape of the gyro and cal- 
culating the integral in equation (2). For a perfectly spherical gyro 
T =0, as is obvious from equation (2). 

To calculate "£ wc need to make severnl assumptions, the validity 
of which we now discuss. 

a) The rotor shape is symmetrical about its spin axis. This 
certainly appears reasonable if the rotor is spinning sufficiently fast, 
since the irregularities would then average out over many rotations of 
the rotori Honeywell [1] has produced an argument which shows that 
this is, in fact, true if one has a sufficiently smooth rotor (see Section 
XI). There is one important exception, however, and that is for a 
rotor with radial mass unbalance perpendicular to the spin axis. In 


r ■ > - — * A A dd I 


this case, the suspension forces cnn have an oscillatory component at 
the spin rate frequency and the radial mass unbalance torque docs not 
average to zero. There is then present a torque component parallel to 
the spin axis which ulters spin speed and can cause run-down. 

b) The rotor shape can be represented by a cosine series of the 

form 


r(lf) * r 0 ♦ 2] n n 008 < n0 ') 
n=2 


(3) 


Almost any shape can be so represented, and the number of terms 
needed is related to the dimensions of the radial irregularities as dis- 
cussed in Section 111. A derivation that does not assume this is given 
in Section V. 

c) The electric field at a given point in the rotor electrode gap 

is given by |£ | = V/d, where d is the rotor electrode distance at that 

point und V is the rotor electrode potential. The justification for this 
is thut tlic rotor electrode gap is so small (1500 gin. typical) that one 
cim treat the situation locally as a parallel plate in which the electric 
field is uniform. Honeywell (1] has made some calculations that bear 
this out. 

d) Electrode edge effects can be neglected. It is difficult to 
estimate how good this assumption is, but it is probably much better 
for hexahedrad electrodes than circular electrodes. 

e) The vaariutlons in rotor electrode gap Ad defined by d = d Q + 
Ad satisfy the conditions 


Ad 

d 

o 


<< 


(4) 


This is a good aissumption. since l lie main contributor to Ad is the 
centering errors which aire of order 15 |iin., giving 


Ad 

d 
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0.01 
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Given all these assumptions, one can write the expression for $ in 
the form 


c r ^ / 

$ _ o o v* « 2 c rl* 2Ad\ dr(e’) ,? , it 

T ■ ^ ‘ J J( " *7/ ~w~ ‘ 8 " ♦ 

o i S| 

- fcos d>*) sin 0' d0' d<t>* , (5) 


where the sum is taken over the six electrodes and 0', <p' are spherical 
coordinates defined about the gyro spin axis. This follows directly from 

2 

equation (2). We have used equation (4) and expanded the 1/d term, 
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set r = r Q in the term dS = r sin O' d0' d<j>', and expressed the surface 

integral as a sum of integrals over electrodes at potentials Vj. The 

expression for (r * n) is justified by noting that it must be tangent to 

a latitude line; i.e.. in the direction given by the unit vector f sin - 

j cos and of magnitude equal to the (small) angle between the radius 
vector and the normal vector to a curve dr/dS* (this is worked out in 
elementary calculus). 

The integral in equation (5) is in a spherical coordinate system 
defined by the rotor spin axis tipped at an arbitrary angle with respect 
to the electrode axes. Since we must calculate the surface integrals in 
a coordinate system defined by the electrode axes, we need to transform 
equation (5) to this coordinate system. If ( 0 Q , $ ) represents the 

orientation of the spin axis in these coordinates, then 

i = cos 0 „ cos *|> F + cos 0 sin 4 > F - sin 0 F 
o ox o T oy oz 

i = - sin F + cos $ F 
J ox o y 

cos «!>' = -r~ -tt frin 0 cos o cos (J> - $ ) - cos 6 sin 0 1 
» sin 0 o T o o 

sin [sin 0 sin (^ - 4> 0 )] 


define the relationship between the primed and unprimed angles. The 
first two equations follow from the fact that the vector in the spin axis 
direction is 
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K = sin 0 cos 6 F 
o ox 


♦ sin 


"o ** 


$, f 

r y 


4 COB 


e o $_ 

o z 


and f is defined to hove no F component , while 1 , f , and £ must be 

Z '■ 

orthonormal. The second two equations follow from expressing (F . F . 

a 4^4 x y 

F ) in terms of (i , j , x) and looking at the expansion of an arbitrary 
z 

vector in the two bases. Substituting these relations into equation ( 5 ), 
we find 



' o r o 






Yl< 0 > 


F y <’'>o 


- •%> * F * 


•o'] 4 


dr(o')/d(o») 


sirTo* 


sin o do d<t> 


( 6 ) 


where u = sin o cos «K |> = sin o sin <f», and y = cos o and similarly for 
v y Q . This is the basic equation wo will use in ull further calcu- 
lations. All that nsods to be done is to take an expression for Ad and 

dr( 0 ')/d 0 ' and culculntc the integrals. Also note that f • K = 0 . The 
important quantity in the expression is dr( o') /d( o’) , which determines 
the magnitude of the torques. This completes the derivation of the 
basic torque equation. The remainder of tho report is devoted to 
evaluating the equation. 
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III. EXPANSION OF THE ROTOR SHAPE IN HARMONICS 


The basic integral in equation (6) has been evaluated at great 
length by Honeywell [1] for hexahedral electrodes. This has been done 
by expanding Ad and drfd’l/de' in a Fourier cosine series and calculating 
the integral separately for each harmonic. The results of this will be 
discussed in Section IV, but here we wish to discuss the convergence of 
such a series. To do this we have takei. the latest roundness measure- 
ments of the best available rotors by Rank-Taylor-Hobson (December 20, 
1976). These measurements are taken over three great circles on the 
rotor; and, since we are assuming azimuthal symmetry, we need only 
half of each great circle. The three curves we have selected are shown 
in Figures 1, 2, and 3, which give Ar as a function of polar angle 6. 

While some improvement may occur due to spin averaging (see Section XI) , 
we believe this is a reasonable place to start in evaluating the harmonics 
for realistic rotors. 

Fig; es 4, 5, and 6, show the Fourier coefficients for the shapes 
given in Figures 1,2, and 3, respectively. These have been obtained 
numerically from the actual assumed shapes. (The first harmonic has no 
significance for these calculations.) These curves show that we do get 
convergence, but it may require up to 20 harmonics depending on the 
scale of the irregularities in the rotor. For example, in Figure 3 more 
harmonics are required because of the sharpness of the bumps. In fact, 
there is an inverse relationship between the number of harmonics required 
and the smallest angular scale on which the irregularities occur. A rotor 
with very sharp bumps will require many harmonics to faithfully approxi- 
mate its shape. The contribution to the integral in equation (6) from 
dr(6')/d@* is proportional to nsin(n6) for the nth harmonic, so we have 
the possibility of increasing contributions to the torque as a result of 
increasingly higher harmonics. In fact, if Reference 1 is studied, it 
will be found that for the higher harmonics the number of terms in each 
torque expression actually increases with harmonic number, os do the 
coefficients of each term. 

These consideration.' illustrate the basic problem with this type of 
calculation. It appears that contributions resulting from higher harmonics 
will be significant when added up, but it becomes extremely tedious to 
calculate equation (6) for the higher harmonics. This is only feasible 
up to about the sixth harmonic and is extremely difficult even up to this 
level, as a perusal of Reference 1 will show. The only way out of this 
dilemma appears to be the method we use in Section V, where we evaluate 
equation (6) numerically for some actual rotor shapes. 

















IV. TORQUES IN TERMS OP HARMONICS 


The method used to evaluate equation (6) is to expand dr(n')/do' 
and Ad in a Pouricr series of the form 


Ar = V, a„ cos(no') 
n 

n=2 


from which dr(0')/d0' is obtained by differentiation, and the part of Ad 
due to rotor nonsphcricity is obtained by Ad = -Ar. The torques due 
to just the dr(0')/d0' term (Ad = 0) are known ns primary torques. We 
have calculated these for circular electrodes up to the fifth harmonic; 
the results are shown in Table 1. The y and z components of T are 
obtained by suitable permutation of u o , y q nnd the P.'s and M.'s [11. 

Honeywell has calculated these torques up to the sixth harmonic for 
hexnhedrnl electrodes. Wc have developed a computer program to 
evaluate our expressions to sec the angular behavior of the torques. 

This is necessary because heretofore the torques have been estimated 
from expressions valid only near the z-nxis, and it is of interest to sec 
how these torques behave for arbitrary orientations. Also, with our 
computer calculation wc can compare the actual magnitudes of the second 
versus fourth and third versus fifth harmonics to see if convergence is 
occurring, figures 7 through 10 show the relative magnitude of the 
torques for the second through the fifth harmonic in arbitrary units for 
typical values of the harmonic coefficients. The second harmonic is small 
because wc have assumed equal preloads. It is clear that the fourth and 
fifth hurmonics produce torques of the same order of magnitude as the 
second and third. A large number of calculations that are not shown 
here bear this out. This is also obvious from the inspection of Tabic 1. 
The expressions in Table 2 arc computed from general expressions for 
the torque as in Table 1 by specializing the expressions to the case 

Y 1, a << 1. l< << 1 and defining 12 = T/Im, where T is the magni- 
o o o 2 

tude of the torque, I = 2/5 mr is the moment of inertia of the rotor and 
w = v/r is the angular velocity of the rotor. The connection between the 
M's and P's in Table 1 and the acceleration f and preload h are derived 
as follows. Prom equation ( 1) we have for circular electrodes that the 
force on the rotor is 


P = 


2 


n sin 



P + M 
x y 


P + M P ) 
y /. z 


( 7 ) 


this determines the components of the acceleration as 


TABLE 1 
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The expression in Table 2 for the aligned third harmonic then agrees with 
that of Table 1 except for the -1/5, which represents the displacement 
of the center of mass for a pear-shaped rotor of this type. Similarly, 
for the misaligned third harmonic, the terms in Tables 1 and 2 agree if 

2 2 

we take @ 0 B 0 and neglect and . This brinps out the fact that 

the expressions in Table 2 ore only good for Y 0 ^ 1 (gimbaled gyro) and 

shows why it is necessary to write r. computer program to evaluate these 
torques for the nongiiubelcd case. This is described in Section V for 
all the harmonics. 

The preload hj is dcllncd here us the acceleration that will drive 

the voltage on one of the ith electrodes to zero and thus lose the 
suspension. (This definition fMTers from that of Reference 1.) For a 
linear suspension system we h. ve V. + V. V. = constant, so from 
oq tuition (7) 1+ *' lo 


mh. 


(V. 2 ) (i. sin 2 u ) 
2 d 2 *° 1 


Using the identity 


(V * ♦ v 2 ) = - 

lv i- v i+ } 2 


< V i* ♦ V>“ + 


( v i+ 2 - V,. 2 ) 2 


< V i* + V 


r 


i? 

2 


1 + 


CV |+ 2 - V,_ 2 ) 2 


T 


io 


* 2 2 
and multiplying by r r / 2 d .we obtuin 

o o o 


p i ! 


m 


2 i! sin o 


h i * 


l L 


if]- 


( 8 ) 


i 
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This is the basic equation connecting Pj with hj and fj. The preloads hj 

are adjusted to be larger than the maximum acceleration expected in a 
given environment. The misaligned second harmonic in Table 2 is seen 
to agree with Table I if is again identified with > and i. is the 

difference in preloads for different axes. 

Comparing these expressions with Reference 2 indicates that there 
are terms present for the aligned and misccntered third harmonic thnt 
are ivu included in Reference 2. Theso two terms put more severe 
restrictions on f and h than have previously been indicated. We post- 
pone generating actual numerical tolerances until Section VI and refer to 
Reference 3, where numerical values for the parameters in Table 2 arc 

* 16 

given which allow the reduction of ii to - 1.6 • 10 rod /sec -• 

0.001 arc scc/yr, the accuracy necessary for a 2 percent measurement 
of the hense T birring effect. 

Tables 3 and 4 contain misccntcring torques for the second ami 
third harmonics in terms of arbitrary electrode angle The correspond 

ing expressions in Table 2 are obtained for the ease 1 as before. 

These expressions illustrate the rapid increase in complexity for higher 
harmonics. 

It is inqiorlanl to realize that the simplified expressions in Table 2 
which have been used in previous drift estimates arc only valid for near 
alignment of the rotor spin with the electrode axis. When this is not 
the case, the expressions in terms of harmonics become extremely tedious 
to evaluate for the higher harmonies, and it is best to resort to numerical 
integration of equation (6). as in Section V. The planned orientation lor 
the final gym experiment is not tlte gim baled case but the case where 
the spin axis is midway between two electrodes. Thus, the expressions 
in Table 2 arc not really relevant for this orientation, and what we 
really need is the case y - t» u \l/F. This is difficult to uvaluatc for 

all the terms in Table I. hut we immediately realize that the alignment 
factor o for tlte second harmonic is no longer applicable. This is a very 
important (mint ami increase i the corrcs|Ntnding drifts by at least two 
orders of magnitude. Tltis |>oint is discussed further in Sections VI 
and Vll. 
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TABLE 3. SECOND HARMONIC MISCENTERING TORQUE (CIRCULAR ELECTRODES) 




TABLK 4. THIRD HARMONIC MISCENTERING TORQUE 
(CIRCULAR ELECTRODES) 
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• 5 

sin o do 




V. EXACT NUMERICAL CALCULATION 


OF TORQUES 


P° s ® il ^ le to evaluate equation (6) for arbitrary rotor soin 

suhitlt!.?" by - d< S n ? the inte K ra,s numerically. To do this we use the 
substitutions m Reference 1 to include all the electrodes; that !?• 


z- electrode 

a -> 

Y -*■ -y 

x+ electrode 

a 

Y. Y -► -a 

x- electrode 

u ■> 

~Y. -*■ -U, y 

y+ electrode 

a 

~ H H v . y -► 

y- electrode 

u ► 

0 . 0 ► - y , y > 


One notices then that permuting u K. and y is equivalent to permutin, 
o’ L o’ and > 0 1,1 expression lor O' and at the same time permutinj 
the three basic integrals given by 


*1 ( V l o’ 

>o> = 

2 ii 

/ ii* 
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r "' . 

J sln 

o 

Odo r„ drop/do'l 

L sin o 1 J 



2 n 

. (, i 

1 sin 
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l 2 (a o* 'o’ 

>o> = 

/ *M' 

o 

0 do [& dr ( l) !2/tio' 
L sin O' 



2/i 

0 

, 1 

J s ’ n 
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! 3 (t V 1 b 1 

V - 

/ <M< 
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u do Tr 

L sin o’ J 

COS O' - UK 
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+ »' ,; o 

+ "o 




a =’sin 0 eos .)> 


= sin o sin 1 j> 
y = eos 0 


25 


Thus, it is possible to express equation (6) in terms of these three 
basic integrals with the arguments permuted in the appropriate fashion. 
The end result of this is as follows (where we have neglected the 
Ad/d term): 


T x - ». |V zl *2 ‘W'.» * V z- '2 

2 d Q { 

* v x? '2 ‘-VW - v x 2 '2 * v y+ 'a 


V v 2 >3 <-vv i; .>" ♦ 'o 1 V W'W 


, 2 ‘3 < VV 'o’ ' V J •, ‘'VW * V *1 ‘-W-..’ 


+ V 


♦ V >1 ( ' V 'l ' 


W'V 1 ) 


r r 

T = -° ° 


y 2 d 


1V J *3 < 'o l o* 'o > V ! 3 
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- V 'l V 3 'l ‘-'o- 'o'-V " V 'l ‘“'o’ 


L \ 

r 2 | 
y+ 2 


‘o» 

* V 

' V J 'i ‘vW 


-> 

>3 < 'V'W * V x 

' l o ) 

•> 

V “ 
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'2 < V V ' o' 1 } 
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T = 


e r 
_o o_ 

2d 2 


b o ' v 2 ; >. <v'w 


■l ^VV^o 1 


+ V x+ ‘3 ( -V’ li o" , o > - V x- '3 ( -V W - V y+ '2 <^o--'*o- l b ) 

+ V 2 ^VVV * •%, '- v ,* >2 <V W 
- V z 2 *2 ‘"VV'V - V x 2 '2 ‘‘V'W * V x 2 '2 ( "V' C V"'o > 


‘ V y t ‘3 ‘‘V'W f V y 2 '3 '"VVV 1 


This is the final exact expression lor the primary torques (with no 
Ad/d Q terms) in terms of an arbitrary rotor shape r(O'), and it does not 

rely on a harmonic expansion. The V. 's can be expressed in terms of f. 

and h. directly from equations (7) and (8) as 



(9) 


We have written a computer progriim (see Appendix) to evaluate 
these integrals numerically for an r(i'’) which we can take directly from 
a roundness chart. The program does the double integrals using an 
8-point Newton Coates quadrature formula with 97 divisions. It has been 
checked with the analytic expressions for the second and third harmonies 
in Table 1 and found to be accurate to four signifie.ant figures. The 
r(O') is Fourier analyzed (0 • 1,1 11 ). anti the Fourier cosine series 

coefficients up to the 20th are used to reconstruct dr(0’)/dt)'. Thus, we 
expect harmonies up to the 20th to be accounted for. 


VI. 


NUMERICAL RESULTS 


A number of computer runs were undertaken to establish the drift 
rates for various values of the parameters involved. The three sample 
rotors were used as previously described. To establish the preload 
dependence, a series of computations were made for a rotor orientution 
of Q q = 45°, <J> o = 0°, where () Q is the angle from the electrode axis and 

<p Q is the azimuthal angle. 

The results are shown in Tables 5 through 8 for various levels of 
acceleration. The h's were taken to differ by 1 percent. The results 
indicate that the drifts are proportional to the preload for the lower 
values of f. The drifts indicate the general order of magnitude to be 
expected for arbitrary rotor orientation. The general level of drift was 

~ 6 

several hundred milliarc scc/yr for 1) ^ 10 g. 


TABLE 5. DRIFT RATES 
(0 o = 45°. -f o - 0°. f •». I0 _, ° g. Ah = 0.01 h) 


h(g's) 

Case 1 

Case II 

Case III 

(milliarc 

sce/yr) 


10 -« 

86.0 

291.0 

202.0 

,0 7 

8.6 

29. 1 

20.4 

io~ 8 

0.89 

3.01 

2.30 


TABLE 6. DRIFT RATES 


<°o 

= 450. * o - 0°. f v. 

10‘ 9 g. Ah = 

0.01 h) 

h (g's) 

Case 1 

Case II 

Case III 


(milliarc 

scc/yr) 


10 -6 

86.0 

291.0 

204.0 

10~ 7 

8.9 

30.0 

23.0 

io“ 8 

2.5 

8.3 

8.5 
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TABLE 7. 

DRIFT RATES 



(0 o 

= 45°, 

♦ 0 =0°, f ' 

- 10 9 g, Ah = 

0.01 

h) 

h(g's) 


Case I 

Case II 


Case III 



(milli arc sec/yr) 



if* 


89.0 

300.0 


210.0 

10" 7 


25.0 

84.0 


85.0 

10- 8 


150.0 

505.0 


370.0 



TABLE 8. 

DRIFT RATES 



<°o 

= 45°, 

<P o = 0°. f -v 10~ 7 g, Ah = 

0.01 

h) 

h (g's) 


Case I 

Case II 


Case III 



(milliar 

c sec/yr) 



io- 6 


254.0 

839.0 


850.0 

O 

1 

^3 


1497.0 

5053.0 


3705.0 


Next the relative effect of the higher harmonics was investigated 
by running the program for each harmonic separately normalized to a 
harmonic coefficient of a^ = 1 iiin. One sees from Table 9 that the 

higher harmonics do not decrease in importance very quickly. Combining 

this with Figures 4, 5. and 6, we see that the a above the second do 

n 

not drop off very quickly, so that the higher harmonics may be important 
and should be included in the calculation. 


TABLE 9. DRIFT RATES (f <», 10 9 g, h 'x. 10" 6 g, a n = 1 uin.) 


Harmonic Number 

Ah = 0.01 h 

Ah = 0 


(niillinrc sec/yr) 


2 

77.0 

0.03 

3 

20.0 

20.0 

4 

19.0 

0.004 

5 

18.0 

18.0 

6 

' 50.0 

0.01 

7 

30.0 

30.0 

8 

16.0 

0.006 

9 

9.7 

9.7 

10 

56.0 

0.01 

11 

11.0 

11.0 

12 

15.0 

0.003 

13 

3.1 

3.1 
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Finally, the factor r, was set equal to zero for the higher harmonics, 
and this greatly reduced the torques for the even harmonics, as shown 
in Table 9. This means that the factor t is generally applicable and 
that the expression for fourth harmonic in Table 1 is probably wrong. 

We have not used it for computational purposes, and we include it in 
Table 1 to illustrate the complexity of the calculations involving higher 
harmonics. 

To summarize, these results indicate that for the rotors considered, 
— 6 

a pretoad h 10 g (determined primarily by roll rate) , and the orienta- 
tion planned for the final experiment (0 O = 45°, <j> Q =0), the drifts are 

of order several hundred milliarc sec/yr, two orders of magnitude larger 
than the experiment goal of 1 milliarc sec/yr. To reach the design goal of 
the experiment one must do one or a combination of three things: 

(1) reduce the preload, (2) align the rotor with the electrode axis (see 
Section VI), or (3) rely on roll averaging of torque. 

While item 2 is a valid solution, it is probably too difficult to 
implement given the increased difficulty of reading out the gyro due to 
the presence of electrodes in the readout loop. Item 3 is also probably 
a valid solution. This is discussed in Section X. 

We notice that for the second harmonic term in Table 1 and a = 

Y q = 1//5 we have 


i2 


X 

Ia> 




1 


2ir sin 0 


l 


(2*7) (a„) . 

( 27 5 T r ~ IT “o Y o a 150 milKarc sec/yr 


for a„ = 2 pin. (due to centrifugal distortion: see Reference 4), h = 

-6 * 

10 g, c = 0. 01 , o> = 200 cps. So just the second harmonic gives torques 
of the magnitude shown in Tables 5 through 10. 

The current state of the art in electrostatic gyros is of order 3 * 
10° milliarc sec/yr in a 1-g environment. Scaling up the results of this 
6 8 

section by 10 , we obtain drift rates of 10 milliarc sec/yr in a 1-g 
environment. Thus, the state-of-the-art rotors considered here will not 
have the 1-g performance of existing electrostatic gyros. This is 
because the current electrostatic gyros produced by Honeywell rely on 
gimbaling, ball balancing, centrifugal distortion compensation, and 
possibly computer modeling to reduce torques. 
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VII. TORQUE ON A GIMBALED GYRO INCLUDING 
ALL THE HARMONICS 


We can evaluate equation (6) without resorting to a Fourier 
expansion for the case y Q = 1, u q = = 0. In this case we have sin6 = 

sin O'; so, neglecting the Ad/d Q term, we obtain 



c r 
o o 


2 d 



f f i» 

S, 


a V 


dr(0) 

de 


d0 d$ 


For the z+ electrode we get 


i = 

z+ 


o o_ 

2d 2 



2 ii 

J d«J> sin 0 sin 4> 
o 


0^ 2n 

- Fy J do J* d<|) sin 0 cos (J> = 0 

o o 


Similarly, for the z- electrode, making the substitution a -*■ -ay-* -y . 
we have 


f = 0 . 

z- 


For the x+ electrode we have a > y , y > -a; so 


c r 2 J. 

' = ° ° V 2 I 1' 

x+ 2 d 2 x+ I x 

r\ 


2/i 0. 


j* d(j> j" sin () di 


( sin 0 sin <i 
J cos 2 0 + 


drCO'l/dQ 1 


2 2 
sin o sin 


2it o. 


- F J* d4> J sin 0 d 


0 / cos 0 dr(6*)/d 
\ / X) s 2 o + si 


2 2 
sin 0 sin <J> 


■) 


31 


O’ * sin / cos O + sin O sin $ 

This follows since the Jacobian of the transformation is sin 0/sin 

sin 6 = / <x + 3 2 . 

Therefore, making these substitutions gives us the radical in the 
denominator. Similarly, so 




IF I. 
x 1 


"yV 




where 


■i 

*2 


2ir «j 

* / d * / sin o do 

o o 

2 it 

- / d<l> / sin 0 do 

o o 


( sin Q sin 4> drCeO/dO* \ 
J cos 2 0 + sin 2 0 sin 2 <J »/ 

(• 


cos Odr(Q t )/dO t 


/ 2 2 2 
cos 0 + sin 0 sin (j) 



For the y+ electrode 


1 11 and 


a -3 . H *► 7 . r ► -a 



y+ 


e r 2 
e o o „ 2 

* V v - 

2 d 


^ y+ 


F x / d* /' dn 0 rtiY ”* 11 dr < f>l > /d11 ' 1 

o o \/ cos 2 o + sin 2 o sin 2 <J> / 


2ii 0j 

+ F J* d<}> J sin 0 d0 j 


( sin d sin 4> drOH/dQ* \ 
J cos 2 e + sin 2 0 sin 2 4> / 


For the y- electrode u > nnd fl ► -y ; so clearly 


$ = -f 

y- y+ 

and 

V M y l V! t Vi' • 

The 1^ integral can easily be shown to be zero by the following argument. 
The <(> part of the integrul cun be written as 

2ir it 

J* f(sin 2 sin <J> di|> = J f(sin 2 <J>) sin 4> d<J> 
o o 

2 ii 

+ f f(sin 2 4>) sin <J> d4> ; 


letting <f> = ii + U in the second integral, then sin 4’ = -sin U and the 
latter integrul becomes 
§ 


TT 

/ f(sin 2 U) (-sin U) ell) 


o 


which exactly cancels the first integral. The final form of the torque 
is then 


* ■ '2 *■"» * ", *x' • 


From Reference 1, page 187, this can be written 


f — I-ft x ftl . 

n sin 0| 


where ft * F is a unit vector in the direction of the spin axis and ft is 
z 

the force on the rotor given by equation (7). This is of the form of a 
mass unbalance torque for the component of unbalance along the spin 
axis. It includes all the harmonics through the integral Ig, which can 

be evaluated exactly for any given rotor. This will lead to a drift rate 
of magnitude 





2 

it sin 2 a* 


m f 

I -V* 


B Lj 

2ir sin Oj 



f 

V 


Using 

f » 10’ 9 g 

V * 2400 cm /sec for w* 200 cps 

r ■ 0.75 in. 0, * 30° 

o 1 


we find 


il * (1.77 x io" 9 ) (I 2 ) 

for Ig in pin. Ig has been evaluated numerically for the three sample 
rotors considered earlier, and the results are 
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Case I 0.983 x to"® in. 

Cose II -0.710 x io ' 6 in. 

Cose 111 -0.995 x uf 8 in. 


Case I ft * 
Case II ft * 
Case III ft * 


10.6 milliarc sec/yr 
7.88 milUnrc sec/yr 
11.0 milliarc scc/yr. 


The drifts are proportional to f and 1/m. So for a perfectly aligned 

-9 

rotor we get drift rates of the order 10 milliarc sec/yr in a 10 g 
environment. 

The large reduction in torques for the gim baled case illustrates 
the great advantage for this orientation. This is, in fact, the orientation 
used for many existing electrostatic gyros and is responsible for their 
current excellent drift performance. These gyros use a rotor with a 
preferred spin axis so the axial muss unbalance torque can be minimized 
by balancing the ball. In fact, all the even harmonics vanish for this 
orientation (1]. As soon as the rotor spin axis is misaligned with the 
electrode aids, one picks up all the even harmonics, and these involve 
the relatively large preloads h. This accounts for the difference in the 
drifts calculated in this section and those calculated in Section VI. We 

can reduce the average f by drug- free control to the 10~ 10 g level, but 

h cannot be so reduced because of the roll of the spacecraft (h ^ 10 8 g 
for 15 min roll period and 10 cm rudius). 




VIII. 8BC0NDARY T0RQUB8 


The secondary torques in by definition the torques that arise 
Ikon the Ad/d 0 term In equation (6). There are four eouroes of second 

ary torques (1): 

a) Rotor aspheridty: 

Ad ■ - Ar(fr) , 


where 


4r ■ r ■ r. 

o 


6* is given by 


y » cos" 1 (oo^ ♦ B8 0 ♦ YY 0 ) . 

so it is a complicated function of the direction oostotes ouB.y and <^.0 O *Y O . 

b) Rotor miscentering: 

Ad « - »* 0 - 6y c - Y* c . 


where . 


<VV*c> 


indicates the vector from the center of the housing to the center of mass 
of the rotor. 

c) Electrode assembly errors: 


^ £ A. 

Ad ■ c -p (a ♦ 6 ♦ 2y) + -p. (a - 8) 
n n 


e 
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where A^ it the translation of the upper hemisphere parallel to the pro- 
jection of F f on the separation plane and A f is the translation perpendic- 
ular to this direction, c * tl for the upper half and c ■ -1 for the tower 
half. 

d) Electrode asphericity: 

Ad ■ R(a,6,y) , 


where R sped flee the deviations from perfect spheridty of the electrodes. 
This requires a model for R. Honeywell (11 has computed the torques 
for the (a) and (b) cases up to the sixth harmonic and in case (c) for 
the second and third harmonics for hexahedral electrodes. Essentially 
no work has been done on cose (U). Extracting numerical values from 
these formulas is difficult and, ns we have seen previously, using only 
the lowest harmonics is mislending. 

The secondary torques could be computed exactly using the method 
described in 8ection V. One defines the three integrals 


1 


1 


( v w 


2tr 

/ d* 

o 


/ sin 
o 


0 dO 



a dr<ew*l 
sin o’ J 


1 


2 


<V e o’V 



V <V p o’V • 


2ti 

/ 


/ sin 0 do 


[ 


Ad . drfo’J/dO’l 
Y tin ft 1 ' J 


where Ad*is a function of .»,3,y and a , 3 . y . For the four previously 

mentioned cases, a complete expression for the torques could then be 
generated by suitable permutation of the direction cosines and the I’s. 

Each case would have 36 terms. Note, however, that, except for case 
(a), upon carrying out the permutations listed in Section V, the l'a 
would not transform into themselves: thus, one would not have Just three 
basic integrals. For cases (b) and (c), one would have six basic integrals 
(Ad proportional to either u.i:, or y); and for case (d), one would have 
twelve basic integrals. This is a formidable computational task that hits 
not been carried out us yet. 
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One important observation about the secondary torques is that 
none of them averages due to spacecraft roll. This is because they con- 
tain a, and y, and this destroys the averaging discussed in Section X 
even for the case of the rotor spin aligned perfectly with the spacecraft 
roll axis. 

In lieu of a complete calculation of secondary torques, we can make 
the following observation. All the secondary torques are of order Ad/d Q 

times the primary torques. So if we take the general level of primary 
torque drift to be of order 100 milliarc sec/yr, which would only require 
a slight rotor improvement, we see that Ad/d Q need only be 1/100 to 

reduce the secondary torques to the milliarc sec/yr level. Since d Q = 

1500 yin. , this requiem a Ad % 10 yin. Therefore, this is the specifica- 
tion we will require on Ad. It is clearly satisfied for case (a). For 
case (b), a centering accuracy of 10 yin. is attainable, but it must be 
remembered that this is an absolute centering accuracy and not the stabil- 
ity of the centering point. The 10 yin. specification also applies to the 
housing centering errors and the electrode sphericity. A more definitive 
statement than this would require an enormous amount of numerical work. 
However, notice that these conclusions agree with the miscentering terms 
in Table 1, which are down by a factor t/d Q from the primary torques 

for the lower harmonics. Also, the calculation in Reference 1 of electrode 
assembly errors bears out that they are down by a factor A/d for the 
lower harmonics. 


There is, however, a qualification to this statement because of the 
fact that the secondary torques, while down by a factor Ad/d Q , may not 

contain the reductions due to symmetry that occur in the primary torques. 

This is illustrated by the miscentered third harmonic term in Table 2, 

which contains no £ term and is of order 15 milliarc sec/yr for a. = 

-6 -2 * 

1 yin., h = 10 g, t/d Q ^ 10 ; and w * 200 cps. The higher harmonics 

will add to this, and the total magnitude of the drift requires the numer- 
ical integration described previously. (The miscentered second harmonic 
gives an insignificant drift for these same parameters and a = 2 yin. 

-9 ® 

and f ^ 10 g.) However, if the miscentering vector remains fixed with 

respect to the electrodes as the spacecraft rolls, then these miscentering 
torques will average to zero because they will cause no change in elec- 
trode voltages as the spacecraft rolls. These considerations indicate 
that a complete numerical integration for secondary torques is probably 
necessary. 


IX, ORBITAL AVERAGING OF GRAVITY GRADIENT FORCES 


The purpose of this section is to compute the average gravity 
gradient acceleration as a function of distance from the center of mass of 
the satellite. The effect of roll averaging will then be computed. 

The gravity gradient acceleration t is given by the equation 


Nvp , 


where V is a matrix of the form 



3 2 V 



and V is the gravitational potential, p is the vector position of the point 
at which the acceleration is cl 'puted with respect to the center of mass. 
Taking V = GM/r, one easily finds that for a rectangular coordinate 
system 


V= gm/ 


3x 2 -r 2 

3xy 

3xz 

3xy 

3y 2 -r 2 

3yz 

3xz 

3yz 

3z 2 -r 


) 


We now consider a circular orbit in the y-z plane with the roll axis in 
the z-direction. This implies the guide star is in the orbit plane, a 
situation that can be realized with a suitable choice of orbit. Then for 
a particular point in the spacecraft specified by p one has 


p_ - const, 
z 

p = p cos Utt 
H y o 

P X = P 0 sin (1,t 
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where u is the roll angular frequency. The orbital position is given by 
y = R cos Qt 
z * R sin fit 
x = 0 

where ft is the orbital angular frequency and R the orbital radius. Sub- 
stituting this into the equation for V and t , one obtains 


GM 

n 5 


- R p o sin wt 

R 2 (3 cos 2 lit - 1) Cp n cos fait) + (3 R 2 cos Qt sin Qt) p 

O 2 

■ i 

(3 R 2 cos fit sin S2t) (p cos wt) + R 2 (3 sin 2 fit - 1) p 

® z 


We can now compute the average values of these quantities, noting that 


< sin (Ot > = . 

10 T 


where T is the averaging time. Similarly, 


^ 2 n4 . . „ _ sin mT . sin (2fi - w) T 

cos fit cos fait > - 2 wT + 4 (2fi - u> ) T 


sin (2fi + t>i) T 
4 (2fi + fa) ) V 


< sin fit cos fit cos ut > = 1 ~4 " 


. 1 - cos (2fi + u) T 
+ 4 (2fi + fa)) t ' 
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< cos fit sin fit > = S - ^Y ' T 


< COS (lit > = 


sin u) T 
uT 


< sin 2 


fit 


> 


1 1 sin 2 fiT 

5 1 S2T 


We see that all components of <t ■ are of order 1/fiT = 2u/# of orbits or 
l/w T - 2 tt/# revolutions and thus average to zero with the very important 
exception of the second term in <f >. Ignoring the term proportional to 
1/fiT, we have 


> 'V' 


GM 


2R 


3 p z 
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For p ^ 10 cm this term is of order 10 g. We see from Section VIII 
z 

that this will make a large contribution to the torque; thus, every effort 
should be made to minimize it. 


o 

Notice that there are <f. > terms in the expression for torque, nnd 
none o:‘ these averages to zero. However, these occur in the combination 

2 -g -6 

<f. >/h.. For gravity gradient forces of order 10 g and h. ^ 10 g, 

* * _ 1 Q 1 

these terms are of order 10 g, the same order of magnitude as the 
<f.> terms. 

To summarize, we have shown that all gravity gradient accelerations 
average to zero under the combined effects of orbital and roll averaging 
with the exception of the component along the roll axis, and this can be 
minimized by placing all the gyros in the plane of the drag- free proof 
mass perpendicular to the roll axis. 


For an elliptical orbit, one can write 


'a . 
X 

V 


- p < — — * sin uit> 

° r 

<f > 

y 

- GM 

„ o < 3y^r_ cog , i)t > + < Y* > p z 

r r 

<f > 
L z J 


V 7 3z 2 -r 2 

( » o < y 5 cos ait > + < 5 > P z J 

P t* 
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where the average is now over an elliptical orbit. One can express y 
and z and thus r as a power scries in the eccentricity e: 


? = cos sit + 
a 


1 

2 


e (cos 2 Sit - 3) 


+ (3 cos 3 fit - 3 cos flt) . . . . 


\ = sin Sit 
a 


+ 


1 

2 


e sin 2 Sit 


2 

+ (9 sin 3 Sit - 15 sin sit) + . . . 


/! 2 
r = / y + z 


a - semi- major axis 


One can convince oneself then fnat 


< 


1 

7 


sin wt > -► 


0 



cos uit > -► 0 


< ^ cos (o t > -* 0, 


just as in the circular orbit case. This is because one will have terms 
of the form 


< cos m (k sit) sin” (t sit) cos wt > , 
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where m, n, k, and t. are integers and these terms average to zero. 
So. after averaging, we obtain 


r<f >i 


•o 

x 



<f > 

= GM 

< > p 

y 


r 5 ** 



„ 2 2 

<f > 


< 3z ; r > p 

. 2 J 


r 5 * J 


But < > is of order e Z /a 3 for an elliptical orbit; so <f > is of 

2 r 3 y 

order e GM p /a , which is small for a nearly circular orbit. Also. 


< 



> is different for a noncircular orbit. 


The corrections are 


r 

order e for small e. So by minimizing p we essentially eliminate the 

z 

the average force in a noncircular orbit also. 


of 


While the z component of gravity gradient acceleration does not 
average to zero, it will not cause substantial torque on the gyro aligned 
with the z-axis. This can be seen by the following argument. The vol- 
tages needed to cancel the z-component of acceleration are such that for 
each point on the rotor the voltage is the same as for a point opposite 
on a line drawn through the spin axis. The torque due to the original 
point thus precisely cancels the torque due to the opposite point. This 
can be seen explicitly for the third and fifth harmonics in Table 1 which 

vanish for a = y =1/ JT , 3 = 0 , M =M,M = 0 , the case appro- 
o 'o ’ ^o ' x z’ y 

priate for the acceleration and spin aligned with the z-axis. The even 

harmonics will still be there, but these are primarily due to asymmetries 

in the preload thnt would be there in any case. The non-zero <f^> for 

an elliptical orbit will cause a torque on the gyro aligned with the z-axis. 
This acceleration is of order 


^ * 



2 

e p 


z 
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and can be reduced to the 10 g level for e ^ . 1. Thus, the non-zero 
average value of the z-component of gravity gradient acceleration will not 
significantly affect the gyro uligned with the z-axis from which the rela- 
tivity information is extracted. ^ 
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X. AVERAGING OF TORQUES DUE TO SPACECRAFT ROLL 


lit the final version of the gyro experiment, the spacecraft will be 
rolled about an axis through the guide star. This is necessary to 
achieve the required readout accuracy, but it is also effective in reducing 
the primary electrical torques on the gyro. In this section, we evaluate 
the efiect of averaging on the torques analytically. 

Consider first the case in which the gyro spin axis is perfectly 
aligned with the roll axis. >Ve know that the torque is perpendicular to 
the roll axis. Then, if everything else is constant, the tip of the torque 
vector will describe a circle, and the average of the torque over one 
revolution will be precisely zero. 

To prove this analytically we note that equation (6) is also valid 
in a coordinate system with the roll axis as the z-axis. Now consider 
an element of electrode surface described by the direction cosines 
a, 3, y. As the spacecraft rolls, these will vrry in such a way that 


a = sin 0" cos ai t 


0 = sin 0 " sin w t 
y = cos 0" 

where 0" remains constant. The rotor spin axis is at the point midway 
between two adjacent electrodes; in the present coordinates this means 



= 0 . 


Hence O' =0” = constant. Thus, if the Vj's remain constant in equation 
(6) then averaging over one roll we have 


<a> = 0 <3> = 0 


Since this holds for each element of electrode area, we have (neglecting 
Ad/d Q terms) 

<f> = 0 
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for one complete revolution. At any given point the torque is reduced 

by a factor 1/n of its static value, so that for 10^ revolutions (%1 yr 
for 15 min roll period) the torques are reduced by four orders of magni- 
tude. This holds only if the V^'s remain constant, which is the case for 

the centrifugal accelerations due to the roll. It is not true for other 
accelerations, but these are taken care of by the drag- free proof mass 
and the gravity gradient acceleration averaging as described in Section 
IX. The averaging also does not hold for the secondary torques (Ad/d Q 

terms), since these terms may contain, a, 0 and y and thus destroy the 
averaging. An important example of this is miscentering torques (see 
Section VIII). 

If the rotor spin axis is not precisely aligned with the roll axis, 
then a V 0, 0 i 0 and we get terms which do not average to zero, 

since yet and yfl terms are present, and also 0 now varies as the 
o o 

spacecraft rolls. These terms arc very difficult to compute, but we sec 
they are of order a Q or p times the primary torques. For example, for 

a 10 arc sec misalignment of the rotor spin axis, a Q and are of order 

5 x 10 so the nonaveraged part of the torques are reduced by this 
factor relative to the static primary torques. 

Thus, we conclude that roll averaging is extremely effective in 
reducing torques provided the spin axis can be kept aligned closely to 
the roll axis. This statement applies only to the primary torques and 
requires that the preloads remain nearly constant during one roll period. 

We can thus sec that the effect of rolling the spacecraft is to elim- 
inate the first h. term in equation (9); therefore, we have run our 

program with this term set equal to zero. The results arc shown in 
Table 10. They indicate that the drift rates can be reduced to the milli- 

-9 

arc sec/yr level for f 10 g for the rotors considered here. We also 
note that the higher harmonics are not very significant. This was also 
the case for the torques before averaging when the relatively large 
second harmonic due to centrifugal distortion dominated the torque 
expression. This term is of order 2 yin. for w = 200 cps and was 
included in the previous calculations of Tables 5 through 8. We conclude 

that milliarc sec/yr drifts arc feasible for f ^ 10 g after averaging due 
to roll is .accounted for. Roll averaging of secondary torques could be 
accounted for also by setting h. = 0 in equation (9) and performing the 

integrations described in Section VIII. 


* 
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XI. SPIN AVERAGING 


Honeywell [1] has given a proof that if the rotor is spinning suf- 
ficiently fast, then one can assume that the rotor is axially symmetric 
with a shape given by 


2ir 

? (0) - ^ f r ( 0,<J>) dip, 
o 


where r (0,<t>) is the actual ball shape in terms of 0 and <p. This is valid 
provided 


271 



where 


A(9,<f>) = r(0,4>) - r(fl). 


Reference 1 refers to this as a "locally smooth" condition. To obtain the 
appropriate input for the torque calculation, one must compute r(8). This 
requires a complete map of the ball in 0 and <p, which is not currently 
available. For the purposes of this report, we have used r(0,<|> o ) for 

0 < 0 < tt and some p instead of r(0). Undoubtedly some averaging 

occurs in the integral of r(0) that will improve the effective smoothness 
of the ball, but it is difficult to estimate how much until a complete map 
of the ball is available. A numerical integration to determine r(0) will 
then be required. 


XII. CONCLUSION 


The main conclusion of this roport is that Newtonian electrical 
torques permit drift rates of order several hundred milliarc sec/yr with 

preloads of 10 g and existing rotors. The preload setting is determined 
by safety considerations, and the roll rate, and correspondingly better 
drift rates are possible with lower preload settings. Roll averaging will 
reduce the primary torques to well below the 1 milliarc sec/yr level. The 
secondary torques are at the 1 milliarc sec/yr level for existing rotors, 
assuming a fractional uniformity of rotor electrode gap of 1/100. The 
gravity gradient forces average out when combined with the spacecraft 
roll except for the component along the roll axis. This component, how- 
ever, does not produce significant torque on the gyro aligned with the 
roll axis from which the relativity information is extracted. The effects 
due to spin averaging of the rotor asphericity may lower the drift rates 
calculated here, but this calculation awaits a complete map of the rotor. 
The overall conclusion is that, from the point of view of electrical torque 
calculations, an accuracy of 1 milliarc sec/yr drift is a feasible goal for 
the gyro experiment. 
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